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We calculate the 0(£)-term of the two-loop massive operator matrix ele- 
ments for twist 2-operators, which contribute to the heavy flavour Wilson 
coefficients in unpolarized deep-inelastic scattering in the asymptotic limit 
S> m?. Our calculation was performed in Mellin space using Mellin- 
Barnes integrals and generalized hypergeometric functions. The 0(e)-term 
contributes in the renormalization at 3-loop order. 

PACS numbers: PACS numbers come here 

1. Introduction 

In the range of small values of the Bj0rken variable x, the contribution of 
heavy flavour corrections to deep-inelastic structure functions is of the order 
of 20-40% and hence has to be known in the QCD analyses of the structure 
functions for high precision extractions of the parton densities and the QCD 
scale Aqcd [!]• In the full kinematic range, a semi-analytic result for the 
heavy flavour Wilson coefficients up to next -to-leading order exists [2] , with 
a fast implementation to Mellin-space given in [3], whereas a fully analytic 
result to O(a^) could be achieved in the limit ^ in [4], Q2 denot- 
ing the virtuality of the exchanged photon and m? the mass of the heavy 
quark. The corresponding Wilson coefficient for at 0(af) was calcu- 
lated in [5]. In this limit, the heavy-flavour contributions can be expressed 
as a convolution of light-flavour Wilson coefficients and massive operator 
matrix elements (OMEs) between light partonic states. The results in [4] 
have been obtained using integration-by-parts techniques. We performed a 
first recalculation of these OMEs in Mellin-space [6,7], using both Mellin- 
Barnes integrals and generalized hypergeometric functions. This shifts the 
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problem of solving complicated integrals of Nielsen-type in [4] , to the cal- 
culation of sums over products of harmonic sums [8, 9] depending on the 
Mellin-parameter N, weighted binomials and Euler Beta-functions. The 
expressions in our result are even on the diagrammatic level considerably 
smaller than the ones obtained in [4], and are seemingly more suitable to 
the problem. In this paper, we show a first step towards the 0(Q;f)-term 
of the heavy-flavour Wilson coefficients, by calculating in dimensional reg- 
ularization the 0(e)-term of the two-loop OMEs, with e = D — 4. 

2. Method 

Our calculation is performed in the asymptotic limit » m^, apply- 
ing the light-cone expansion, where, as the massless renormalization group 
equation (RGE) gives a splitting of the deep-inelastic structure functions 
F2/L into a convolution of perturbatively calculable Wilson coefficients and 
non-perturbative parton distribution functions, the massive RGE allows to 
write the heavy flavour contribution to the twist-2 Wilson coefficients as a 
convolution of light -flavour Wilson coefficients and massive operator matrix 
elements [4]: 

massive OMEs light Wilson coefficients 

These OMEs are universal objects, calculable via the corresponding flavour 
singlet, pure-singlet and non-singlet operators between partonic states, de- 
termining the non-power contributions in m^/Q^. The process dependence 
is then solely given by the massless light Wilson coefficients [10]. 
The OMEs contain ultraviolet and collinear divergences, the former be- 
ing removed through renormalization, the latter absorbed into the parton 
distribution functions. To two-loop order, the renormalized gluonic OME 
reads: 




with similar expressions for the quarkonic contributions. Here, P^j are 
the fcth-loop splitting functions, /?o is the lowest order expansion coefficient 

9 (k) 

of the /3-function, and /x the renormalization and factorization scale, a^j 
and a^^^ are the 0{£^) resp. 0(£:)-terms in the expansion of the OME. As 
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Fig. 1. Two example diagrams contributing to the gluonic OME, with aU fermion 
lines massive and external momentum = Q. 

a first step towards a 0(af) calculation, one needs each of these quantities 
to one additional order in e, since they then enter the constant term of the 
OME by multiplying the corresponding splitting functions. The 0(e)-term 

(2) _(2) 

of the OME A- , a- , is a new result presented here and the main topic of 
this calculation. 

3. Calculation 

The diagrams can be grouped into two sets: one-loop in one-loop insertions 
and generic two-loop diagrams. They are calculated using FORM [11] and 
MAPLE programs. Figure [1] shows some diagrams contributing to the gluonic 
OME. 

The rules for operator insertion are, e.g., given in [12]. The calculation 
is done on the one hand by the use of Mellin-Barnes integrals to produce 
numeric results. These results serve as a check for the analytic results, ob- 
tained by expressing the diagrams as generalized hypergeometric functions 
which are first expanded in e and then summed up to the desired order. 
The application of Mellin-Barnes integrals for scalar diagrams in our frame- 
work has already been explained in some detail before, cf. e.g. [6] (see 
also [13]). The idea is to express in a loop-by-loop manner the sub-diagram 
of a full diagram into a Mellin-Barnes representation and to combine this 
with the remaining part. For full diagrams with a numerator structure, one 
can make heavily use of the fact that the light -like vector A occurring in 
the numerators obeys = 0. This reduces the integrals to be calculated 
to a smaller set. After finding a suitable Mellin-Barnes integral representa- 
tion, we use the mathematica package MB [14] to numerically calculate the 



Table 1. Mellin moments N= 2 and N= 6 for diagram e. 



Diagram N 






1 


e 


2 

£ 


e 2 


8.88889 


-11.2593 


9.82824 


-12.8921 


2.39145 


6 


2.93878 


-4.24257 


3.39094 


-4.3892 


0.826978 
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Mellin-Barnes integrals for fixed values of Mellin-A^, up to a given order in 
e. As an example, Table [1] shows up to some generic multiplicative factors 
the results for the full diagram e of Figure [TJ 

By closing the contour and applying the theorem of residues, it is in 
principle possible to even obtain analytic results from the Mellin-Barnes 
representation [6]. However, the way of hyper geometric functions turned 
out to be more appropriate for the calculation of analytic results. In this 
case, one first introduces Feynman parameters and does the two momentum 
integrations. As an example, the scalar version of diagram e of Figure [1] with 
all propagators to the power one, evaluates to: 



(Ap)^^-ir(l-e) 



iV(iV + l)(47r)4+^(m2)i-£ 



dz dw ■ 



w 



(z + W 



wz 



W 



N+1 _ (1 _ ^)JV+1 



On rewriting this Feynman parameter integral into a p-Fg-function, one ob- 
tains a product of 3i<2-functions and the Euler Beta-function 
B{a,b) ■.= T{a)T{b)/T{a + b): 



. i=2 



' N{N + 

B{e/2 + l,l-e/2)B{l,-e/2) 3F2 
- B{e/2 + 1,1- e/2)B{l, N + l- e/2) 3F2 



l-e,l,l + e/2 ' 
2,1 -e/2 



- B{e/2 + 1,1- e/2)B{N + 2, -e/2) 3F2 



l-e,l,l + e/2 
2, N + 2 -e/2 ' 

1 -e,iV + 2,l + e/2 
2,iV + 2 - e/2 



c : = 



(47r)4(m2) 



One then expands this expression up to the desired order in e, in this ob- 
taining finite and infinite sums over harmonic sums and Beta-functions: 

. g.fl_5^^5.(A^ + ,)_B(/V+M)1 
1) , s s s \ 



iV(iV + i)^ 



It is the next step of summing up the various sums over harmonic sums 
and more complicated expressions, which constitutes the most difficult part 
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of the calculation. These sums could be solved among other things using 
their integral representations, where a certain amount of more complicated 
sums could be calculated using the mathematica package SIGMA [15,16]. For 
integral e, one obtains up to 0(e): 

c \ SI{N) + ^S2{N) Sl{N) + ZSi{N)S2{N) + SS^{N) \ 
" A^(iV + l)| 2 12 

In a similar manner, it was possible to calculate all diagrams contributing to 
the calculation of the 0(£)-term of the two-loop unpolarized OMEs. Here 
algebraic relations between harmonic sums were used [17]. 

4. Results 

The 0(e) contributions to the mass-renormalized unpolarized OMEs for the 
singlet, pure-singlet, and non-singlet cases read: 



7t(2) - t r [2 (iV^ + iV + 2)(3JV^+3jV + 2) , Pi „ 

«Qs ^«^^|3 N\N + l)\N + 2) ^'^iV3(iV+l)3(Ar + 2)''' 

JV4 _ 5iv3 _ 32jv2 _ i8jV - 4 ^ + 2 (^^^ 

iV2(iV + l)2(iV + 2) + N{N+l){N + 2) l^^^^'i'i - 



+ - 



— S&.ift + 3yS4 — -S3S1 — —S2 — S2S1 — -Si + 2(^282 — 2(^2Si — -CsSi) 

OA D o ^ 



iV^ - 3iV - 2 ^ 2 3Af + 2 ^3 



m{N +1){N + 2) ' 3m{N + 2)' 
2 3N* + 48jV3 + 43JV^ - 22JV - 8 3iV + 2 Sj_ 

"^3 iV2(iV+ l)2(Af + 2) N^{N + 2) m''^ 

N^ + N'^- 8N'^ -5N^ -3N -2 

2N^ - 2N* - IIN^ - 19iV^ - 44JV - 12 P2 

iV2(iV + l)3(7V + 2) '"^ iVS(iV+l)5(iV + 2) 



f N'^ + N + 2 / 
+^^^-^1 N(N + l)(N + 2) (16^-2.1.1 - 4^2.1.1 - 85-3.1 - 8S_2,2 - 4S3.1 - 



+954 - leS-a.ift + y ftSs + 4/3" 51 - 8/3' S2 + ^Sl - 8/3' 5? + SS?^ + ifi^ 



-ySiCs - 2S2C2 - 2S?C2 - 4/3'C2 - yd) 



■i{N^-N-A) / „ 2 A/"^ + 8iV2 + llAf + 2 , 



- (AT + 1Y{N + 2)2 V-^^-^'i + ^" - ^'^'^O ~ 3 iV(iV + 1Y{N + 2)2 
AT^ + 2Ai'3 + 7iV^ + 22Af + 20 , 3JV^ - 12A^^ - 27iV - 2 
^ (iV+l)3(iV + 2)3 ^ Ar(iV + l)2(7V + 2)2 ^ ' 
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16 N'' + lOA''^ + 9iV^ +3N^ + 7N + 6 N'^ + N - 1 

(iV- l)iV2(Ar+ l)2(Af + 2)2 ^ ^ (A/-+ l)2(iV + 2)2^''^ 

2 9Af^ - lOiV* - lliV^ + 68Af^ + 24Af + 16 P3 

~3 (iV- l)iV2(iV + l)2(iV + 2)2 (N -1)N^{N + 1)^(N + 2)^ 

2P4 P5 ^2 

{N - 1)N^N + l f{N + 2)2 N{N + l)3(iV + 2)3 ' 

2P6 2P7 



iV(Af + l)4(iV + 2)4 (N -1)N^{N +1)^N + 2y 



' = TfC^|-54 + -S2C2 - ^SiCa - y53 - y5iC2 + 2"g^(^";^)"C3 
112 ^ SiV* + eAf^' + 47N^ + 20iV - 12 . 656 „ Pg 

+ T:Z-'52 H 1 n Ar9 / ar , i \9 ^2 ^-Jl + 



4^ 20 „ 20 „3Ar2 + 3Ar + 2 



27 18iV2(Ar + l)2 81 648Af4(7V + l) 

-ps,(2) / , (5iV^ + 7iV^+4iV + 4)(jV^+5jV + 2) ^^^ 

^Q, - TfC^ < -2 (^_i)^3(;v+l)3(iV + 2)2 1^52 + (2) 



(iV2+7V + 2)2(3S3 + C3) 



3(iV-l)iV2(iV+l)2(Ar + 2) (TV - 1)7V^(7V + l)^(iV + 2)4 

with the polynomials Pi given by 

Pi =3iV'' + 30W'^ + ISiV* - 64Ar^ - 56N^ - 20iV - 8 , 

P2 =8iV^'' + 24iV^ - UN^ - 128iV^ - 195N'^ - IWN^ - 23iV* - 27Af^ - 45iV^ 

- 24iV - 4 , 

P3 =A^'' + 21iV* + S5N'' + lOSiV*^ + 42Ar'^ + 290Af'' + 600iV^ + 456iV^ + 256Af + 64 
P4 =(iV^ + 3iV^ + 12iV + 4)(iV'^ - iV* + 5iV^ + iV + 2) , 
P5 =iV'* + 6Ar^ + 7iV'* + 4N^ + 18iV^ + 16A'' - 8 , 

P(i =2N^ + 22Af^ + 117iV® + 386iV" + 759Ar'' + 810iV^ + S96N^ + 72iV + 32 , 
P7 =4iV^^ + SOiV^" + 267Af" + 765iV" + 1183iV" + 682Af"' - 826iV^ - 1858Af* 

- 1116Af^ + 457Af'^ + 1500iV^ + 2268Ar^ + 2400A''^ + 1392iV^ + 448iV + 64 
Pg =1551A/'** + 6204Ar^ + 15338iV'^ + 17868iV'^ + 8319iV'' + 944iV^ + 528Af^ 

- 144Ar - 432 

Pg =5iV" + 62iV"^ + 252iV^ + 374iV** - AOON^ + 38N'' - 473iV^ 

- 682iV* - 904Ar-' - 592iV^ - 208iV - 32 . 

Here, 5'i,...,j = Si^,,,j{N), (3 = (3{N + 1), and Ci = C(*) is Riemann's Zeta- 
function. In all results, there is an overall factor Sla1{m? / fj^Y and an 
overall factor (1 + (— 1)^)/2 in all singlet and pure-singlet cases. The 
above Mellin-space expressions can be converted to x-space using analytic 
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continuation in N and Mellin inversion [18]. 

5. Comparison 

The results in [4], which are up to constant order in e, involved 48 basic 
functions, cf. [7]. Our results for the constant part of the OMEs gave raise 
to only six harmonic sums, where five of them can be obtained from Si{N) 
through differentiation after analytic continuation and using algebraic rela- 
tions. This leads to an amount of only two basic functions [19]. To order 
0(e), we encounter the following 14 harmonic sums: 

{Si, S2, S3, S4, S-2, S-3, S-4}, 52,1, 'S'_2,l, S-S^i, 52,1,1, 'S'_2,l,l, 

S-2,2, 53,1. 

We can again group the first seven functions into the same class. Addition- 
ally, one finds that the function 5_2,2 depends on 5_2,i and 5_3,i, and the 
harmonic sum 52,1 depends on 53,1, which leaves us to order e with only six 
basic harmonic sums, as also observed for a large variety of other two-loop 
processes, cf. [20]. 

6. Conclusion 

We have calculated the 0(e)-term of the unpolarized massive two-loop 
OMEs, contributing to the heavy-flavour Wilson coefficients in the asymp- 
totic limit ^ m^, as a first step towards the 0(Qf)-term of these Wilson 
coefficients. The calculation was done in Mellin space, where numeric results 
were obtained by the use of Mellin-Barnes integrals, whereas analytic results 
were calculated using generalized hyper geometric functions. After applying 
algebraic relations, the analytic result for the 0(e)-term is expressible in 
only six basic harmonic sums. 
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